Long range rapidity correlations and the ridge in A+A collisions 
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Abstract 

We discuss results for n-gluon correlations that form the basis of the Glasma flux tube picture of 
early times in heavy ion collisions. Our formalism is valid to all orders in perturbation theory 
at leading logarithmic accuracy in x and includes both QCD bremsstrahlung and the many body 
screening and recombination effects that are important at large parton densities. Long range 
rapidity correlations, as seen in the near-side ridge in heavy ion collisions, are a chronometer of 
these early time strong color field dynamics. They also contain information on how radial flow 
develops in heavy ion collisions. 

Recently, high energy factorization theorems were derived [ lj for inclusive multi-gluon pro- 
duction in a rapidity interval Ay < l/a s in A+A collisions. This work is based on a field theory 
formalism for particle production from strong time dependent sources J2]. Our result can be 
expressed very compactly as HI 
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In the Color Glass Condensate formalism [3], the Z's are gauge invariant weight functionals that 
describe the distribution of color sources in each nucleus at the rapidity of interest and are ob- 
tained by evolving the JIMWLK equations [4] from an initial rapidity close to the beam rapidity. 
The JIMWLK equations contain information about multi-parton correlations at all rapidities in 
the nuclear wavefunctions. In the large N c limit, for large nuclei, the Z's can be obtained from 
the simpler mean field BK [5] equation. 

The leading order single particle distributions in Eq. (HJ are given by 
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For each configuration of sources p \ 2 of each nucleus, one solves classical Yang-Mills equa- 
tions to compute the gauge fields A&\p\,pj\ in the forward light cone Jg]. Equation (13, when 
substituted in Equation (HJ, gives, to all orders in perturbation theory and to leading logarith- 
mic accuracy in x, the n-gluon inclusive distribution in high energy A+A collisions. Eq. (HJ is 
valid when all the produced particles are measured within a rapidity interval < l/a s from each 
other; Z is evaluated for rapidity y\ j» • ■ • » y n ta y. The expression in Eq. ([T} suggests that in 
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a single event-corresponding to a particular configuration of color sources-the leading contri- 
bution is from n tagged gluons that are produced independently. The correlations in n-particle 
gluon emission are generated by event averaging over color sources. Because the range of color 
correlations in the transverse plane is of order 1 / Q$, where Qs is the saturation scale, our re- 
sult suggests an intuitive picture of coherent multiparticle production arising from fluctuations in 
particle production from color flux tubes of size 1 /2s- 

These color flux tubes, called Glasma flux tubes J3l, are approximately boost invariant in 
rapidity because of the underlying boost invariance of the classical fields. The Glasma flux 
tubes also carry topological charge [8], which may result in observable metastable CP-violating 
domains [9]. Albeit eqs. ([]} and (0 describe initial particle production from the Glasma flux 
tubes, they do not describe the growth of instabilities or subsequent final state interactions of the 
produced gluons which become important for times t > l/2s- Issues related to the "unstable 
Glasma" and its decay have been discussed by us elsewhere 11 1 CHI . 

At large transverse momenta {p ± > Qs), analytical expressions for the single particle distri- 
butions in eq. © are available [11]. These were used to compute 2-gluon [ 12], 3-gluon 1 1311 and 
n-gluon correlations [14] in the MV model tJ. One obtains 1 14] 
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where k = ( n (N 2 - l)Q^S ± /2n, („ is a non-perturbative coefficient that can be determined from 
numerical simulations, and Q\S ±_ ~ N part . The expression in Eq. (01 corresponds to a negative 
binomial distribution which is widely seen in hadronic scattering data. The variance in these 
distributions can be expressed as <x 2 = n 2 - h 2 — h + h 2 /k. Both a 2 and k are proportional to 



Npart, as seen in the RHIC data 11511 . This remarkably simple geometrical structure of multi- 
particle correlations was derived for p ± > Qs ; however, there is now evidence from numerical 
lattice simulations of two particle correlations that this structure may be robust at low p ± as well. 

The ridge [16] in the Glasma flux tube model is a combination of initial state long range 
rapidity correlations and final state radial flow which collimates the angular distribution of cor- 
related pairs. From causality, the latest time at which gluon rapidity correlations could set in 
is given by the relation t < tj exp(-0.5Ay), where Ay is the rapidity separation between pairs 
and Tf is the freeze-out time. The persistence of the ridge to large Ay suggests that rapidity 
correlations are formed at early times; indeed they reflect correlations in the nuclear wavefunc- 
tions as suggested by eq. (HJ. The angular correlation can be seen already from simple radial 



boost models where the boost parameter is fitted to the single particle spectra 112LI13I1 . A more 



sophisticated treatment gives a good fit to RHIC data on the centrality dependence and energy 



dependence of the amplitude of the ridge and of its azimuthal width II 1711 . Computations II 1311 . 
are also consistent with observed features of the data on three particle correlations. 

The blast wave model of radial "Hubble" flow needs to be improved to account for all global 
features of the ridge data. An example of this is the "away side" ridge, which likely is due 
to global energy momentum conservation in each event. A first hydrodynamical calculation 
incorporating flux tube like structures in the initial conditions reproduces all such features of the 
data 01811 ; however, only "triggered" ridge events have been considered thus far. With a realistic 
hydro model, ridge like correlations, when compared to data, could provide further insight into 
the development of radial flow and the persistence of long range correlations at late times lfl9ll . 

We now turn to genuine long range correlations. As stated, eq. (HJ is valid for Ay < l/a s , 
which corresponds approximately to 3-5 units of rapidity at RHIC and the LHC respectively. The 
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only energy dependence in this expression comes in from the energy dependence of the weight 
functionals Z (equivalently, Qs), while the dependence on Ay is flat. This is an approximation 
because gluon radiation occurs in the rapidity interval between tagged gluons. Further, gluons 
are emitted in the presence of colored sources; the inclusive multiparticle spectrum at large Ay 
is therefore sensitive to non-trivial multi-parton correlations in the nuclear wavefunctions. This 
problem has been solved B20I1 : the solution is represented diagrammatically in Fig. ([T]). The 
modification to eq. (HJ is that one now has in addition Green functions describing the evolution 
from gluon y p to y q and vice-versa. These Green functions obey the JIMWLK equation. As the 
rapidities of the two gluons approach each other, one recovers eq. (Q}. Phenomenological work 
exploring the consequences of this novel formalism will be reported shortly. 




Figure I: Left: Gluon ladder diagram demonstrating building blocks in QCD evolution of the inclusive 2-gluon spectrum 
at large rapidity separations. The merging of conributions of several sources is typical of JIMWLK evolution. The crosses 
on rungs of the gluon ladders denote the positions in rapidity of the two tagged gluons. The evolution from nucleus 2 
at the bottom of the figure is inferred. Right: Relative magnitude of leading order near-side Long Range Correlations 
("Glasma")and near-side Short Range ("pQCD") Correlations. Color Charge density of projectile 1 (p\ = 0(g~')) is 
held fixed while that of target 2 (p2) is varied. See text for details. 

We turn finally to discuss the relation of our work on Glasma correlations to correlations 
from perturbative QCD ("pQCD"). Examples of respective leading order diagrams and their 
power counting are shown in the right figure of Fig. (Q]). The Feynman diagram on top depicts 
two particle correlations from disconnected graphs in the Glasma. These generate the near- 
side long range rapidity correlations (LRC). The bottom figure depicts the contribution to two 
particle correlations from 2-to-2 Feynman graphs, analogous to conventional pQCD two particle 
correlations. These generate the near-side short range rapidity correlations (SRC). The middle 
graph is the contribution to 2 particle correlations from the interference between the two. When 
the parton densities in both nuclei are large (central collisions/high energies/low momenta), the 
cross-section for LRC dominates that of the SRC; as suggested by the positions of the filled dots, 
the disconnected graph contribution is larger than the interference graph by 1 /a s and is larger 
than the pure SRC graph (at bottom) by as much as an order of magnitude (1/a^). In addition, 
one has a kinematic suppression of the latter cross-section. This is because "pQCD" correlations 
are SRC for good reason-they die off rapidly with Ay. In contrast, as we have discussed here, the 
disconnected Glasma graphs are a natural QCD mechanism to generate long range correlations. 
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SRC will begin to compete in magnitude with LRC as parton densities are decreased (peripheral 
collisions/pA/high momenta). In pA collisions, as schematically represented by the filled dots in 
the left hand corner of the plot, the two mechanisms are parametrically of the same order, even 
though one expects the relative kinematic suppression of pQCD to persist. In pp collisions, SRC 
dominate, except perhaps at very low momenta. These statements can be quantified for a unified 
("the long and the short") QCD perspective on multi-particle correlations. Work on extending 
our formalism (and phenomenological consequences thereoff) to pA collisions is in progress. 
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